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To understand the pressure-induced changes in the electronic structure and the electron-phonon 
interaction in yttrium, we have studied hexagonal close-packed (hep) yttrium, stable at ambient 
pressure and double hexagonal close-packed (dhep) yttrium, stable up to around 44 GPa, using 
density-functional-based methods. Our results show that as one goes from hep yttrium to dhep 
yttrium, there is (i) a substantial charge-transfer from s — > d with extensive modifications of the d- 
band and a sizable reduction in the density of states at the Fermi energy, (ii) a substantial stiffening 
of phonon modes with the electron-phonon coupling covering the entire frequency range, (iii) an 
increase in the electron-phonon coupling constant A from 0.55 to 1.24, leading to a change in the 
superconducting transition temperature T c from 0.3 K to 15.3 K for ji* — 0.2. 



I. INTRODUCTION 

The change in the superconducting transition temper- 
ature T c in yttrium from 6 mK at ambient pressure^ to 
around 20 K at 115 GPa 2 , concurrent with the change 
in its crystal structure 3,4 , underscores the intricate rela- 
tionship between pressure-induced s — > d charge trans- 
fer, dynamical response of the lattice and its coupling 
to Fermi electrons. The experiment on yttrium fol- 
lows the recent observation of superconductivity in com- 
pressed lithium at a T c of 20 K— bringing into focus 
the pressure-induced changes in the electron-phonon in- 
teractions in elemental metals, especially the alkali and 
the early-transition metals£aSa£. 

Experimentally, yttrium follows the equilibrium struc- 
ture sequence hcp^Sm-type^dhcp^trigonal as the 
pressure is increased from ambient to 50 GPa at room 
temperature 3 -. It is expected that yttrium will stabilize 
in a body-centered cubic structure at pressures above 280 
GPa 10 . On the other hand, The superconducting tran- 
sition temperature T c changes from 6 mK at ambient 
pressure 1 - to around 1.3 K at 11 GPa 11 , and then to 17 
K at 89 GPa 2 . A further increase in pressure to 115 GPa 
results in a T c = 20 K. Surprisingly, it seems possible 
that for pressures greater than 115 GPa, yttrium could 
superconduct with a T c higher than that of 20 K. 

In an effort to understand the pressure-induced 
changes in the superconducting properties of yttrium vis- 
a-vis the changes in its stable crystal structure, we have 
studied, using density-functional based methods, (i) the 
electronic structure, (ii) the phonon density of states, (iii) 
the electron-phonon interaction, and (iv) the solutions of 
the isotropic Eliashberg gap equation for hep and dhep 
yttrium at ambient pressure and at around 50 GPa. 

We have calculated the electronic structure of hep and 
dhep yttrium using full-potential, linear muffin-tin or- 
bital (LMTO) metho d 12 ! 13 as well as the plane-wave 
pseudopotential method using Quantum-ESPRESSO 
package^ 4 -. Using the linear response code based on the 
LMTO method, we have calculated the phonon DOS, 
F(u), and the Eliashberg function, a 2 F(u;), for hep and 
dhep yttrium. Subsequently, we have numerically solved 
the isotropic Eliashberg gap equatio n 15 ! 16 ! 17 for a range 



of fi* to obtain the corresponding superconducting transi- 
tion temperature T c . After the completion of the present 
work, we have come across the work of Yin et alt^, which 
uses the hypothetical fee structure to describe the effects 
of pressure in yttrium. 

Based on our calculations, described below, we find 
that as one goes from hep yttrium to dhep yttrium, there 
is (i) a substantial charge-transfer from s — > d with ex- 
tensive modifications of the d-band, (ii) a stiffening of 
phonon modes with the electron-phonon coupling over 
the entire frequency range, (iii) an increase in electron- 
phonon coupling constant A from 0.55 to 1.24, leading to 
a change in the superconducting transition temperature 
T c from 0.3 K to 15.3 K for /i* = 0.2, consistent with 
experiment. 



II. THEORETICAL APPROACH AND 
COMPUTATIONAL DETAILS 

The density-functional theory provides a reliable 
framework for implementing from first-principles the 
Migdal-Eliashberg approach for calculating the supercon- 
ducting properties of metals. Such an approach has been 
implemented using the full-potential LMTO metho d 12 ! 13 . 

In general, the phonon density of states is given by 



0, 



where (j qry is the qv-th phonon mode associated with the 
atomic displacements. Using the electron-phonon matrix 
elements ^k+qj' kj> which can be interpreted as the scat- 
tering of electron in state \kj > to state |k + qj' > due 
to perturbation arising out of the phonon mode a; qi/ , we 
can calculate the phonon line width 7 qi/ 
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where N(ef) is the electronic density of states at the 
Fermi energy. Now, we can combine the electronic 
density of states, phonon spectrum and the electron- 
phonon matrix elements to obtain the Eliashberg func- 
tion a 2 F(co) defined as 
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Finally, the Eliashberg functions can be used to solve the 
isotropic gap equatio n 15 ! 16 ! 17 

\LJ n f \<u c 

A(i(j n ) = ^2 fn>S(n,n)A(iu n >) 

n' 

to obtain the superconducting properties such as the su- 
perconducting transition temperature T c . The function 
S{n,n r ) used in the gap equation is defined by 

S(n, n) = A(n - ri) - // - 5 nn > ^ s n s n » A(n - n") 

n" 

and f n = l/|2n + 1| with s n representing the sign of c^ n . 
The electron-phonon coupling A(V) is given by 

Before describing our results, we provide some of the 
computational details of the present calculations. 

The self-consistent electronic structures of hep and 
dhep yttrium in P6s/mmc crystal structure, having 
two and four yttrium atoms in the respective primi- 
tive cells, were calculated using the pwscf code of the 
Quantum-ESPRESSO package and the full-potential, 
LMTO method. The results obtained with the pwscf 
code are used to study the equation of state of yttrium. 
All other properties of yttrium, described below, have 
been obtained with the full-potential, LMTO method and 
linear response code based on the LMTO method. 

The plane- wave pseudopotential method with Vander- 
bilt's ultrasoft pseudopotentials were used to calculate 
the electronic structure of hep and dhep yttrium with 
259 k-points in the irreducible Brillouin zones. The ul- 
trasoft pseudopotential used also included the 4s and 4p 
states with nonlinear core correction. The kinetic en- 
ergy cutoff for the wavefunction was taken to be 35 Ry, 
while for the charge density and the potential the kinetic 
energy cutoff were 200 Ry. For exchange-correlation po- 
tential we used the local-density approximation as given 
by Perdew et al 19 . These calculations were used to study 
the equation of state of yttrium. 

The charge self-consistent, full-potential, linear muffin- 
tin orbital calculations for hep and dhep yttrium were 
carried out with the local-density approximation for 
exchange-correlation of Perdew et al^ with 2^-energy 
panels and 610 and 549 k-points respectively, in the cor- 
responding irreducible wedges of the hexagonal Brillouin 
zones. The 4s state of yttrium was treated as a semicore 
state, while the 4p state was treated as a valence state. 
The basis set used consisted of s, p, and d orbit als at the 
yttrium site, and the potential and the wave function 
were expanded up to l m ax — 6. The muffin-tin radii for 



yttrium at zero pressure and at 50 GPa were taken to be 
3.3 and 2.8 atomic units, respectively. These calculations 
were used to study the s-, p- and d-resolved densities of 
states, the band-structure along the high symmetry di- 
rections in the hexagonal Brillouin zones and the Fermi 
surfaces of hep and dhep yttrium. The Fermi surfaces of 
hep and dhep yttrium were constructed using XCrySDen 
program^ with eigenvalues calculated on a 36 x 36 x 24 
and 36 x 36 x 12 grid in the reciprocal space respectively. 

The phonon density of states and the Eliashberg func- 
tion of hep and dhep yttrium were calculated using the 
linear response code based on the full-potential, LMTO 
method. The dynamical matrices and the electron- 
phonon matrix elements of hep and dhep yttrium were 
calculated on a 6 x 6 x 4 and 6x6x2 grids respectively, 
resulting in 21 and 14 irreducible q-points. The Brillouin 
zone integrations during linear response calculations for 
hep and dhep yttrium were carried out using a 12x12x8 
and 12x12x4 grids of k-points. The Fermi surface sam- 
pling for the evaluation of the electron-phonon matrix el- 
ements for the two structures were done using 36 x 36 x 24 
and 36 x 36 x 12 grids. 



III. RESULTS AND DISCUSSION 

In this section we describe the results of our calcu- 
lations of the electronic structure, the linear response 
and the solutions of the isotropic Eliashberg gap equa- 
tion for hep and dhep yttrium. Our results are described 
in terms of (i) equation of state, (ii) electronic structure, 
(hi) electron-phonon interaction, and (iv) superconduct- 
ing transition temperature obtained from the solutions 
of the isotropic Eliashberg gap equation. 



A. Equation of State 

The calculated equations of state of yttrium in hep 
and dhep phases are shown in Fig. 1. These equations of 
state are obtained by fitting the calculated total energy 
versus volume points to the third-order Birch-Murnaghan 
equation of stated. The value of c/a was optimized by 
iteratively minimizing the total energy as a function of 
volume. The lattice constants a and c at zero and 50 
GPa for both hep and dhep phases, as obtained through 
the above procedure, are given in Table I. However, the 
lattice constant a for the high-pressure dhep phase of 
yttrium corresponds to the compressed volume equal to 
0.56^0, with the ideal c/a = 3.25 and where vq is the 
experimentally determined equilibrium volume of the hep 
phase. The lattice constants as given in Table I have 
been used in all the calculations reported below. We 
also note that in our discussion, we will emphasize the 
stable structures, viz. hep at zero pressure and dhep at 
50 GPa pressure, more than the unstable structures at 
these pressures. 
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Figure 1: The equation of state (top panel) of hep (solid line) 
and dhep (dotted line) yttrium, the difference in energy AE 
(middle panel, dashed line) and enthalpy AH (bottom panel, 
dashed line) of dhep yttrium with respect to hep yttrium as 
a function of volume per atom. 



Table I: The lattice constants of hep and dhep yttrium used 
in the calculations. The experimental values 3 are given in the 
parentheses. 



structure 


pressure (GPa) 


a (a.u.) 


c/a 


hep 





6.775 (6.895) 


1.56 (1.571) 


hep 


50 


5.604 


1.65 


dhep 





6.811 


3.20 


dhep 


50 


5.619 


3.25 



The calculated lattice constants for hep yttrium at zero 
pressure, as given in Table I, are somewhat smaller than 
the experimental values, as expected in local-density ap- 
proximation. We find the bulk modulus of hep yttrium 
at equilibrium to be 39.4 GPa, which compares well with 
the experimental value of 41 GPa 3 . In Fig. 1, we also 
show the energy and enthalpy of dhep yttrium with re- 
spect to hep yttrium as a function of volume. At zero 
pressure the dhep phase is approximately 3 mRy higher 
in energy than the hep phase. However, with increase 
in pressure both energy and enthalpy of the dhep phase 
becomes lower than the hep phase, in agreement with 
experiment. 



B. Electronic Structure 

We show in Fig. 2, the calculated band-structure 
along high symmetry directions in the hexagonal Bril- 
louin zones of hep and dhep yttrium at both zero and at 
50 GPa pressure. The band corresponding to 4p state is 
much lower in energy, and it is not shown in the figure. 
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Figure 2: The band-structure along high symmetry directions 
in the hexagonal Brillouin zones of (a) hep yttrium at ambi- 
ent pressure, (b) dhep yttrium at 50 GPa, (c) hep yttrium 
at 50 GPa and (d) dhep yttrium at ambient pressure. The 
horizontal line, passing through the energy zero, indicates the 
Fermi energy. 
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We find that in the case of hep yttrium at zero pressure 
only two bands, the 3rd and the 4th bands in Fig. 2 (a), 
cross the Fermi energy. Based on the band-structure, 
we expect the states close to the zone boundary such 
as M, K, H, and L symmetry points to contribute to 
the electron-phonon coupling. For dhep yttrium with 4- 
atom basis, we find that there are four bands crossing the 
Fermi energy at 50 GPa, as shown in Fig. 2 (b). There 
are states around symmetry points M, K and H that are 
close to Fermi energy, therefore we expect these states 
to play an important role in electron-phonon coupling in 
compressed dhep yttrium. Due to charge transfer from 
s — > d the bottom of the s-band moves slightly up in 
dhep yttrium as can be seen by comparing Figs. 2 (a) 
and 2 (b). 

In Figs. 2 (c)-(d), we show the band-structure of hep 
phase at high pressure and dhep phase at ambient pres- 
sure, respectively. We find that the increase in pressure 
has changed the band-structure of the hep phase, espe- 
cially along T-A and around symmetry points K and H, in 
such a way that it looks similar to the high pressure dhep 
phase. In particular, the flat portion along K-H in the 
hep phase, as shown in Fig. 2 (a), has been completely 
modified at high pressure, as shown in Fig. 2 (c). 

A more comprehensive picture of s — > d charge transfer 
can be obtained by comparing the densities of states of 
hep and dhep yttrium. In Fig. 3, we show the total 
and the /-resolved densities of states of hep and dhep 
yttrium. By comparing the total densities of states in the 
two structures, we find that many more states are created 
between —0.26 Ry and —0.16 Ry in the compressed dhep 
yttrium to accommodate the pushed out s-electrons as 
well as the ^-electrons pushed further inside from the 
flattened portion of the c?-band near the Fermi energy. 
We also find that in the energy interval — 0.2 Ry to —0.05 
Ry, some p-states are also created in the compressed dhep 
yttrium. 

The density of states around the Fermi energy plays an 
important role in determining the electron-phonon inter- 
action in metals. In the present case, we find that the 
total density of states of 22.8 st/Ry atom at the Fermi 
energy in hep yttrium reduces to 13.4 st/Ry atom in dhep 
yttrium. The reduction in s and d states are from 0.71 
and 20.9 st/Ry atom in hep yttrium to 0.6 and 11.7 st/Ry 
atom in dhep yttrium, respectively. Usually, an increase 
in the density of states at the Fermi energy strengthens 
the electron-phonon coupling. However, in spite of a de- 
crease in the density of states at the Fermi energy as yt- 
trium is compressed, there is an increase in the strength 
of the electron-phonon coupling. 

The changes in the densities of states of the unsta- 
ble phases of yttrium, namely hep at high pressure and 
dhep at ambient pressure, are shown in Fig. 3 (b). The 
pressure-induced changes in going from ambient pressure 
dhep phase to high pressure hep phase are qualitatively 
similar to that of hep to dhep as described above. For 
example, the total density of states of 21.1 st/Ry atom at 
the Fermi energy in dhep yttrium, at ambient pressure, 
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Figure 3: The s, d and the total densities of states of (a) 
hep yttrium at ambient pressure (solid line) and dhep yttrium 
at 50 GPa (dashed line), and (b) hep yttrium at 50 GPa (solid 
line) and dhep yttrium at ambient pressure (dashed line). The 
vertical, dashed line indicates the Fermi energy. 



reduces to 13.1 st/Ry atom in hep yttrium at high pres- 
sure. The changes in s and d states are from 0.58 and 
19.3 st/Ry atom in dhep yttrium to 0.6 and 11.5 st/Ry 
atom in hep yttrium at high pressure, respectively. 

The topology of the Fermi surface is useful in under- 
standing many of the reciprocal-space-based properties 
of the solids. In particular, flat portions of the Fermi 
surface can give rise to peaks in the response function, 
leading to a substantial increase in the response function. 
In the case of compressed lithium, it has been shown that 
the flat portions of the corresponding Fermi surface, re- 
sulting in Fermi surface nesting, substantially enhance 
the electron-phonon interaction^^. 

In the present case, the calculated Fermi surfaces of 
hep and dhep yttrium are shown in Fig. 4. In hep yt- 
trium, the 3rd band gives rise to a cylindrical tube along 
T-A, connected to a pancake-like structure with hexago- 
nal symmetry on either end, while the fourth band gives 
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Figure 4: The Fermi surface of (a) hep yttrium at ambient 
pressure, (b) dhep yttrium at 50 GPa, (c) hep yttrium at 50 
GPa, and (d) dhep yttrium at ambient pressure. 



rise to a cylindrical tube along T-A with a bigger radius. 
Due to pressure, the Fermi surface of yttrium is substan- 
tially altered, as can be seen in Fig. 4 where we also 
show the Fermi surface of dhep yttrium. The most ob- 
vious thing to note is the reduction along /^-direction of 
the Fermi surface of dhep yttrium by half. It is due to 
the doubling of the unit cell along the ^-direction as one 
goes from hep to dhep structure. The cylindrical tube 
has separated from the pancake-like surface and, in ad- 
dition, we find some flat portions on the M-K-H-L plane 
with webbing 22 like features. 

As shown in Figs. 4 (c)-(d), the change in pressure does 
not substantially alter the shape of the Fermi surfaces 
of either hep or dhep phase. However, to be able to 
compare the Fermi surfaces of hep and dhep phases at 
ambient pressure, we have calculated the Fermi surface 
of hep phase after doubling the length of its primitive cell 
along c— direction. A band-by-band comparison of the 
two Fermi surfaces of hep and dhep phases at ambient 
pressure is shown in Fig. 5. The smaller values of a and 
c/a in hep phase compared to dhep phase are responsible 
for bigger Brillouin zones corresponding to the hep phase. 
The differences in the Fermi surfaces of the two phases 
can then be attributed to the way the four basis atoms 
are distributed in the respective unit cells. 



C. Electron-Phonon Interaction 

The primary objective of the present study is to un- 
derstand the pressure-induced increase in the supercon- 
ducting transition temperature T c of yttrium. Within 
the conventional BCS theory of superconductivity, a 
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Figure 5: The Fermi surface of hep (left column) and dhep 
(right column) yttrium, both at ambient pressure, corre- 
sponding to (a-b) band number 1, (c-d) band number 2, (e-f) 
band number 3, and (g-h) band number 4. 



change in T c arises due to changes in the phonon spec- 
trum and/or the electron-phonon spectral function (also 
known as the Eliashberg function) . Usually, with increas- 
ing pressure the lattice hardens resulting in an increase 
in the phonon frequencies. At the same time, it is also 
possible for some phonon modes to soften and, in some 
cases, it may lead to lattice instability. 

In Fig. 6, we show the phonon density of states 
F{uo) of stable hep and dhep yttrium, calculated using 
the density-functional theory as described earlier. The 
phonon density of states of hep yttrium, as given in Fig. 
6, is somewhat harder than the results of Heid et a&. 
The maximum phonon frequency is about 29 meV with 
a large peak centered at 26 meV. We find that for q = 0, 
the maximum frequency has a contribution from the two 
basis atoms moving along the £-axis in opposite direc- 
tions. By comparing the F(u) of dhep yttrium with that 
of stable hep phase, we clearly see the hardening of the 
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Figure 6: The phonon density of states F(u) (top panels in (a) 
and (b)) and the Eliashberg function a 2 F(oo) (bottom panels 
in (a) and (b)) of (a) hep yttrium at ambient pressure (solid 
line) and dhep yttrium at 50 GPa (dashed line), and (b) hep 
yttrium at 50 GPa (solid line) and dhep yttrium at ambient 
pressure (dashed line). 



lattice as reflected in the increase of the maximum fre- 
quency to 41.8 meV with a broad peak at around 36.9 
meV. We find that for q = 0, the maximum frequency 
has contribution from the four basis atoms, separated 
into two groups of equivalent atoms, both groups moving 
along the £-axis in opposite directions. 

To see the pressure-induced changes in the electron- 
phonon interaction in yttrium, we have calculated the 
Eliashberg function a 2 F(uo) of hep and dhep yttrium as 
shown in Fig. 6. A comparison of a 2 F(uo) of hep and 
dhep yttrium shows that the electron-phonon coupling 
is weaker in hep yttrium than in dhep yttrium. In par- 
ticular, the average electron-phonon coupling constant A 
is equal to 0.55 for hep yttrium and 1.24 for dhep yt- 
trium. The pressure-induced increase in A occurs due to 
enhanced coupling between 12 — 25 meV and due to the 
coupling of those phonons that arise due to the hard- 
ening of the lattice. Thus, we find that the Eliashberg 




Figure 7: The partial A q along high symmetry directions 
in the hexagonal Brillouin zones of (a) hep yttrium at am- 
bient pressure (filled circles connected with solid line) and 
dhep yttrium at 50 GPa (filled squares connected with dashed 
line), and (b) hep yttrium at 50 GPa (filled circles connected 
with solid line), and dhep yttrium at ambient pressure (filled 
squares connected with dashed line). The solid and the dot- 
ted lines connecting the calculated points are only a guide to 
the eye. 



function a 2 F(u) extends up to 41.8 meV in dhep yttrium. 
The higher frequency couplings are more effective in rais- 
ing the superconducting transition temperature than the 
lower frequency couplings. 

The role played by the crystal structure in determining 
F(uo) and a 2 F(uo) can be seen from Fig. 6 (b), where we 
show F(u) and a 2 F(uo) for the unstable hep and dhep 
yttrium phases. Not surprisingly, now it is the ambient 
pressure dhep phase of yttrium which has a lower maxi- 
mum frequency of 27.5 meV with the high pressure hep 
phase having a higher maximum frequency of 43.6 meV. 
Consequently, the electron-phonon coupling is weaker in 
the dhep phase than in the hep phase, leading to A = 0.47 
for dhep yttrium and A = 1.2 for hep yttrium. 

To gain further insight into the nature of electron- 
phonon coupling, we show in Fig. 7 the q-resolved partial 
A q along symmetry directions in the corresponding Bril- 
louin zones for stable and unstable hep and dhep phases 
of yttrium. We find that as we go from ambient pres- 
sure hep phase to high pressure dhep phase, the electron- 
phonon coupling enhances substantially around symme- 
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Table II: The calculated average electron-phonon coupling 
constant A, the root mean square frequency uo r ms and the su- 
perconducting transition temperature T c obtained with /i* = 
0.2, for yttrium as a function of structure and pressure. The 
T c values given in the parentheses correspond to /x* = 0.12. 



yttrium 


pressure (GPa) 


A 


cj rms (K) 


UK) 


hep 





0.55 


193.5 


0.3 (2.8) 


hep 


50 


1.21 


221.7 


15.9 (19.8) 


dhep 





0.47 


209.0 


< 0.001(1.4) 


dhep 


50 


1.24 


214.2 


15.3 (19.0) 



try points A and M, as can be seen from Fig. 7 (a). The 
contribution to A from symmetry point A is A q = 8.1, 
which arises from the lowest eight modes involving mo- 
tion of all the four atoms in the x — y plane. Similarly, at 
the symmetry point M, A q = 4.1, with the lowest mode 
involving all the four atoms contributing around 1.6. 

A comparison of A q for unstable dhep and hep phases 
with that of the stable hep and dhep phases shows signif- 
icant contributions from symmetry point K, in addition 
to symmetry points A and M. The value of A q = 8.6 
at symmetry point A arises from the doubly degenerate, 
bond-stretching mode of the two basis atoms in the x — y 
plane, each making a contribution of 3.8 to A q . 

We like to point out that the integrated values of A q for 
the unstable phases may be similar to the stable phases 
but the details are very different. We have also seen that 
a substantial contribution to a 2 F(uj) and consequently 
to A comes from a small region in reciprocal space, us- 
ing a uniform grid to sample the space requires a dense 
grid of points. Thus, our calculated values of a 2 F(oo) 
and consequently, A can be improved upon by including 
more q-points in the linear response calculations but it 
is unlikely to change the main conclusions of the present 



work. 



D. The Superconducting Transition Temperature 

The possibility of superconductivity in hep and dhep 
yttrium within the present approach can be checked by 
solving numerically the isotropic gap equatio n 16 ! 17 using 
the calculated Eliashberg function a 2 F(uo). The results 
of such a calculation for hep and dhep yttrium for two 
values of /i* are shown in Table II. At ambient pressure, 
we find hep yttrium to be superconducting with T c = 2.8 
K for fi* = 0.12 and T c = 0.3 K for /i* = 0.20. In the high 
pressure dhep phase, the T c increases to 19.0 K for /i* = 
0.12 and 15.3 K for /i* = 0.20, respectively. In Table II 
we have also listed the electron-phonon coupling constant 
A, and the various averages of phonon frequencies for hep 
and dhep yttrium at both ambient and 50 GPa pressure. 
IV. CONCLUSIONS 

To understand the pressure-induced changes in elec- 
tronic structure and electron-phonon interaction in yt- 
trium, we have studied hexagonal close-packed (hep) yt- 
trium, stable at ambient pressure and double hexago- 
nal close-packed (dhep) yttrium, stable up to around 44 
GPa, using density-functional-based methods. Our re- 
sults show that as one goes from hep yttrium to dhep 
yttrium, there is (i) a substantial charge-transfer from 
s — > d with extensive modifications of the c?-band, (ii) 
a stiffening of phonon modes, (iii) an electron-phonon 
coupling over the entire extended frequency range with 
electron-phonon coupling constant A changing from 0.55 
to 1.24 and (iv) a change in the superconducting transi- 
tion temperature T c from 0.3 K to 15.3 K for /i* = 0.2, 
consistent with experiment. 
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